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Abstract 



We derive a "generic" inhomogeneous "bridge" solution for a cosmological model in 
the presence of a real self-interacting scalar field. This solution connects a Kasner- 
like regime to an inflationary stage of evolution and therefore provides a dynamical 
mechanism for the quasi-isotropization of the universe. In the framework of a stan- 
dard Arnowitt-Deser-Misner Hamiltonian formulation of the dynamics and by adopting 
Misner-Chitre-like variables, we integrate the Einstein-Hamilton- Jacobi equation corre- 
sponding to a "generic" inhomogeneous cosmological model whose evolution is influenced 
by the coupling with a bosonic field, expected to be responsible for a spontaneous sym- 
metry breaking configuration. The dependence of the detailed evolution of the universe 
on the initial conditions is then appropriately characterized. 



1 Introduction 



As is well known [H (21 (see also the general solution of the Einstein equations 

near a cosmological singularity exhibits an oscillatory stochastic behavior. This feature 
of the very early universe is in striking contrast with the universe as described by the 
well-tested theory of the standard cosmological model [Zj, which is based on the highly 
symmetric Friedmann- Robertson- Walker geometry. However, the experimental evidence 
for the homogeneous and isotropic character of our actual universe concerns relatively 
late stages of evolution. Indeed the good agreement of the light element nucleosynthesis 
prediction with the observed abundances implies that the standard cosmological model 
is surely valid after 10^^-10^^ sec from the big bang, but says nothing about the very 
early dynamics before this time. 

In this respect, by observing that the Friedmann- Robertson- Walker metric is an 
unstable solution of the Einstein equations when regarded as running backwards in time 
[Hj, then from the existence of structures in the universe jH] like galaxies and clusters 
of galaxies, we may infer (even in the presence of an inflationary scenario) that such 
symmetric geometry cannot continue all the way up to the initial singularity. In fact, 
the dumpiness of the universe indicates the necessity for very early perturbations of 
its homogeneity and isotropy, which unavoidably "explode" when approaching the big 
ban^ 

The instability of the FRW metric, when regarded backward in time, means that there 
exists some moment t^, before which the evolution of the Universe was to be described 
by a "generic" inhomogeneous model, or , by other words, when the BKL (Belinskii - 
Khalatnikov - Lifshitz) picture |2] holds. In general this moment represents a free 
parameter which depends on initial conditions and, in particular, on specific properties 
of matter. In vacuum case the applicability of the BKL picture was shown [Hj to be 
described by the inequality Lh <^ where L/^ ~ t is the horizon size and Lj„ is the 
characteristic scale of inhomogeneitv (or, to be more precise, the mean geometrical value 
of all leading inhomogeneity scales □) 

Thus, in the vacuum case the moment t^, = tin when the Mixmaster phase (i.e. the 
oscillatory regime) ends corresponds to the situation Lh ~ Lm, which can be roughly 
considered as a boundary of the BKL approximation (in the early stages only the stiff 
matter influences the evolution of the metric). The reversibility of the Einstein equations 

^ Actually, this represents a rather weak point since the universe could start out in a completely 
homogeneous and isotropic state, while primeval perturbations which necessary for the origin of 
structures might appear later during the inflationary era. 

^We recall that inhomogeneity of the metric is described by a set of distinct length scales which 
correspond to different terms in the spatial scalar curvature. The anisotropic evolution leads to 
the increase of some curvature terms (i.e., the corresponding scales decrease with respect to the 
horizon size) which causes the transition between Kasner regimes as described in Ref. 2 . From the 
qualitative standpoint, the transition of a Kasner epoch occurs when the horizon size matches the 
smallest of the scales ^ and therefore the above inequality is violated at least for one of the 
scales. However, the duration of such violations are small compared with the duration of Kasner 
epochs and for the average geometrical value of the scales, this inequality still holds. Anyway the 
greater the inequality, the better the BKL picture works. 
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means that the BKL picture works in both directions of time (which is exphcitly reahzed 
in the so-called billiard representation The fact that tj„ is a free parameter means 

that near the singularity at a moment to, the initial conditions can always be chosen in 
such a way that tin and there is a defined period (which depends on the degree of 
inhomogeneity of the metric) when the BKL picture still works even with the increase 
of time. 

We note that there exists another important moment when the matter "switches 
on" (i.e., starts to influence the evolution of metric). This moment represents another 
free parameter specified by initial conditions for matter. Therefore, the actual moment 
when the Mixmaster evolution (BKL picture) breaks down depends on the relation 
between this two parameters tm and 

In the present paper we consider a restricted region of initial conditions in which 
the inequality tm <^ tm holds, which means that the matter starts to dominate deep 
inside the BKL regime. From the physical standpoint this restriction means that the 
transition from the one BKL regime to another takes place when locally, i.e., on the 
scale of causal connection, the universe still looks like a homogeneous model, for the 
inequality Lh ^ t ^ Un is fulfilled, while global properties are described by a general 
inhomogeneous metric. Therefore, it is natural to expect that some results obtained for 
homogeneous models can also be applied there. In particular, if tm corresponds to the 
beginning of an inflationary period, then the inequality Lh <C tm remains valid from 
the BKL era through the inflationary period H and therefore we can use the results of 
Ref. [Ej. 

The chaotic nature of the evolution (both temporally and spatially) implies that the 
geometry of the very early universe should be described by a stationary statistical distri- 
bution [2] (see also [TTj -[El). Indeed in this context we may speak about geometry only 
in an average sense; it turns out that mean values of all geometrical quantities (lengths, 
scalar products, etc.) during the oscillatory regime are unstable (higher moments have 
the same order of magnitude as the average values) and therefore near the singularity the 
universe does not possess a stable background. We remark that the same situation holds 
in the quantum evolution of the inhomogeneous Mixmaster fSJ universe, although 
in the quantum case the statistical distribution has a different (but somehow related) 
nature. 

These considerations first pose the problem of the origin of a stable background and 
second how this background could arise out of this chaos compatible with the notion of 
isotropy (on the basis of any acceptable early history of our actual universe). However, 
the strong correlation between the appearance of a stable background and its isotropic 
character is a key feature of the very early cosmology. Either on a quantum level or 
on a classical one the isotropic component of the metric tensor (i.e. the volume of the 
universe) is a monotonic function of the time variable (which may actually be taken 
as the time coordinate itself) and therefore does not contain any physical degrees of 
freedom, which are instead entirely contained in the anisotropic components. In other 

■^We recall that during the inflation the inhomogeneity scales increase more rapidly than does the 
horizon size. 
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words a stable background metric can only appear when the anisotropy of the universe 
is sufficiently suppressed jHl HE] • 

In vacuum inhomogeneous models this problem was considered ffist in J7j, which 
outlines how a classical background can arise from the Planckian epoch of the universe 
essentially when the oscillatory regime is over, i.e. at the moment when the characteristic 
scale of inhomogeneity Lin matches the horizon size Lh- Analogously, in the presence of 
matter described by vector fields, the background appears when the horizon size reaches 
the minimal scale between the one related to the vector field and the characteristic 
scale of inhomogeneity ^H] (for a discussion of chaos in superstring cosmology relative 
to Einstein-dilaton-p-form fields see ^H] and [20] )• It is important that in both cases 
the anisotropy of the universe decays (i.e. it becomes smaller and smaller) during the 
natural Mixmaster-like evolution from the initial singularity (it may be worth noting 
the analysis presented in [21] on quiescent cosmological singularities). 

In this paper we consider the origin of a background space when a real self-interacting 
scalar field is present in the universe. As we shall see in this case the appearance of a 
background depends on the initial conditions (to be assigned on a nonsingular space- 
like hypersurface) ; in the configuration (phase) space there are two regions of initial 
conditions for which the evolution behaves in qualitatively different ways. 

The ffist region corresponds to the case when the potential term of the scalar field 
becomes a dominating term before the end of the Mixmaster evolution [Lh ~ Lin) 
(e.g., such a region can be characterized by the inequality Lc <C Lin, Lc denoting the 
Compton length associated with the scalar field, so that the horizon size matches ffist the 
Compton scale). In this case the scalar field through its energy completely governs the 
quasi-isotropization process (i.e. the process which gives origin to a stable background). 
The appropriate region of initial conditions contains a subregion which corresponds to 
an inflationary-like evolution of the universe. The second region corresponds to the case 
when the scalar field potential remains small and, from a qualitative point of view, the 
origin of a stable background occurs in the same way as in vacuum models. 

Below we will consider the ffist region only, having in mind the idea that a classi- 
cal quasi-isotropic universe may emerge, up to suitable initial conditions, from general 
cosmological dynamics, essentially by virtue of an inflationary expansion due to the po- 
tential term of the real scalar field. Indeed the main result of this paper is to show the 
existence of a set of initial conditions of a nonzero measure, corresponding to which the 
anisotropy of the universe decays exponentially during an inflationary phase (in homo- 
geneous models the inflationary phase and the isotropization of the Universe has been 
considered in Ref. ^D])- Thus the analysis of the "generic" cosmological solution shows 
how the inflation phenomenon is the "bridge" between the chaoticity near the big bang 
(indeed in the presence of a real scalar field the Mixmaster contains only a finite number 
of oscillations [221 12S]) and the phenomenology of the standard cosmological model. 

In section 2 we develop the standard Arnowitt-Deser-Misner (ADM) Hamiltonian 
formulation [2ll|2nj which is at the foundation of our derivation in section 3 of a "generic" 
solution of the Einstein-Hamilton-Jacobi equation in presence of a real self-interacting 
scalar field. Such a solution interpolates between a Kasner-like regime and an inflationary 
scenario and is to be regarded as the main result of this paper (for a related discussion 
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of the Bianchi I model in the path integral formalism see P^). In Section 4 we provide a 
reformulation of the system dynamics in terms of Misner-Chitre-like variables, in order 
to give the most appropriate framework for the analysis presented in Section 5 and 
devoted to emphasize the modification induced in the details of the universe evolution by 
assigning different initial conditions to the dynamical quantities involved in the problem. 
Thus this analysis defines the range of existence for the solution we have obtained. 



2 Hamiltonian Formulation of the Dynamics 

Let us start by fixing the dynamical framework for our investigation of the "generic" 
inhomogeneous dynamics. We first observe that the dynamical regime we find will be 
regarded as "generic" in the sense that it possesses the number of physically arbitrary 
functions of the spatial coordinates (i.e. real degrees of freedom of the physical system) 
required to specify a generic Cauchy problem on a nonsingular spatial hypersurface 
(having in mind one which is arbitrarily close to the big bang) for a generic [n + 1)- 
dimensional space-time, containing a real self-interacting scalar field, the number of 
physically independent degrees of freedom is n{n — 1). Indeed this number is (n + 
l){n — 2)/2 for the gravitational field, i.e. the number of independent polarizations of 
a gravitational wave, plus 1 for the real scalar field, but both these fields satisfy second 
order equations. 

The line element of a generic (n-l-l)-dimensional space-time (for the sake of generality 
we will consider the most general case and the results for our actual universe will follow 
immediately by setting n = 3) admits the following standard (ADM) representation: 

ds"^ = N^dt^ - gapidx'^ + N"dt){dx^ + N^dt). (1) 

Then the Einstein-Hilbert action takes the form (in what follows we use units of the 
Planck length) 

1 = j d^xdt + - NH^ - N^h}^ , (2) 

where the super-Hamiltonian if ° and the super-momentum Ha are respectively 

= ^ - iK? + g{W-R)Y (3) 

= -2V/37rf + n<^9„0, (4) 

where W{(j)) = \ {g°'^da4>dji(l) + ^(0)} and R is the spatial scalar curvature constructed 
from the spatial metric gap- 

A fundamental step in our investigation consists of rewriting the above (ADM) for- 
mulation in a form which is useful for our purposes, but which retains its degree of 
generality. It turns out to be convenient to use the so-called generalized Kasner-like 
parameterization of the dynamical variables in terms of n logarithmic scale variables qa 
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and n spatial frame covector fields dual to a spatial frame consisting of n vectors L" 
(having inverse component matrices). The n- dimensional metric components and their 
conjugate momenta are represented in terms of these variables as follows 

g^p = J2 m , = 5^ p^L^t^ , (5) 

a a 

where the covectors £° can contain only n{n — 1) arbitrary functions of the spatial 
coordinates. By definition the Kasner vectors are eigenvectors for both the momentum 
matrix tt^ and the metric Qais and therefore in the case of a generic point in spacetime, 
such a decomposition is unique (n (n — 1) arbitrary functions contained in Kasner vectors 
i'^ and 2n functions and pa replace the n{n + 1) functions contained in ga/s and vr^g). 
This general form of the metric is the most suitable for treating the Kasner-like regime. 

A further refinement of the parameterization can be made by separating the different 
types of contributions to the matrix i'^ as follows (e.g., see Ref. |2j): 

ia = U%Sl U\eSOin), S\ = 5\ + R\ (6) 

with R'^ct denoting a triangle matrix (/?"« = if a < a) and, therefore, it contains 
only n{n — l)/2 arbitrary functions of coordinates, while the rest n{n — l)/2 arbitrary 
functions are included in the rotation matrix U^. By substituting ^ and (jHI) into 0, 
we rewrite the action functional in the form (any repeated index is to be regarded as 
summed) 

^ = / ^^'^^ + + ^^ft - ~ N^H-)crxdt, (7) 

where = '^^^PbL^Ujl and the Hamiltonian and the momentum constraints take the 
form 

E^^-;7^(E^'^)' + ^nJ + f/|, (8) 

a a J 

Ha = -^dp {V9T^S\) + padaq"" + T^daRl+Tl^daCt>. (9) 

We note that due to the property UijdtUik= —UikdfUij the time derivative of the matrix 
(/"'b drops out from the expression (|7j) and besides, only n (n — 1) /2 components of the 
matrix T" (i.e., with a > a) should be considered as independent functions which are 
canonically conjugate to the triangle matrix variables -R"q. 

In the super-Hamiltonian constraint the quantities R'^ and T" contribute only 
to the spatial curvature in the term U = g {W — R) and for the case of n = 3 the func- 
tions R'^ are connected purely with transformations of the coordinate system and may 
be removed by solving the super-momentum constraints if" = [3 (which expresses 
independent components of T" as functions of pa, q"", H^, and 0). In the multidimen- 
sional case, however, the functions R^ contain dynamical functions as well which 
cannot be removed by coordinate transformations. However, in what follows we shall 
use model representations for the potential term U in which the dependence of U on 
i?" and can be neglected (e.g., the generalized Kasner model (GKM) corresponds to 
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the case where we neglect the spatial curvature, the inhomogeneous Mixmaster model 
(IMM) corresponds to the case where we replace the spatial curvature term with a set 
of infinite potential walls). Therefore in these models the super-Hamiltonian will not 
depend explicitly on i?^ and T" and these functions will have a passive character and 
can be considered separately. Indeed, in the case of the GKM or IMM the evolution 
of these functions is completely governed by the supermomentum constraint which 
can be used to express n independent functions among and T" via the rest passive 
[R^ and T") and dynamical (pa, g", 11^, 0) functions. We note that in the GKM all the 
passive functions represent merely constants of the motion (e.g., see Ref. [S]), while in 
the IMM the oscillatory evolution is accompanied by a rotation of the Kasner vectors 
which is completely determined by the momentum constraint (e.g., see Ref. [211221 )• 

After having shown how the above formal decomposition of the metric variables into 
scale functions and "reference" vectors acquires a precise dynamical meaning in the 
above action, we must make a key distinction among the scale functions themselves by 
extracting their isotropic component from the anisotropic ones. This is accomplished by 
further refining the metric parametrization by introducing coordinates on the space of 
scale variables which are quasi-orthonormal with respect to the DeWitt minisuperspace 
metric following Misner |1] 

= A^^(5^ +a /?" = [n(n - l)/2]-^'^(t) , (10) 

where j = l,...,n— 1, and the suitably chosen constant matrix A"- obeys the conditions 

5^A%=0, ^A%A\ = n(n-l)5,,, (11) 

a a 

for example 



_ / n(n-l) fl« - / 1' 3>a, (.^^ 

'~^lijTTj^^~^^ ^ "10, j<a. ^^^> 

Since g = exp(na), we see that a corresponds to the volume or isotropy degree of 
freedom, while (3^ {j = 1,2,. ..,n — 1) describe the anisotropy of the model. When 
these variables are sufficiently suppressed in the sense that they asymptotically approach 
constant values, we may speak of quasi-isotropization of the model. 

Then the action expressed in these variables formally resembles the action of a rela- 
tivistic particle moving in a potential (here the index r runs from 1 to n) 

where the potential term U = n{n — 1)U may be viewed as a "mass term" for the 
particles which depends on dynamical variables (on the position in the phase space) and 
is not everywhere positive. 
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It is important to emphasize that since our goal is to consider simple models (GKM 
and IMM) , for the sake of simplicity we removed from the above action the passive func- 
tions (the terms T^dR^/ dt) as well as the super-momentum terms which as explained 
above can be considered separately. 



3 Construction of the Inhomogeneous "Bridge" Model 



In this section we derive a "generic" inhomogeneous solution connecting the Kasner- 
like behavior (to be regarded as one of the Kasner epochs in the oscillatory regime) 
with an inflationary regime |2ZI- It is worth noting that, although our analysis is done 
on a purely classical level and in the absence of ultrarelativistic matter, nevertheless 
from the qualitative point of view it has a predictive character even in a more general 
context. Indeed on the one hand we may expect that during the last Kasner epoch of 
the oscillatory regime the so-called quantum potential plays, in the Hamilton- Jacobi 
equation, the role of a small correction to the classical potential (which can be modelled 
by a set of infinite potential walls, e.g. see the next section). From the other hand, 
either during the Kasner regime in the presence of a scalar field or during the inflationary 
phase, the contribution of the ultrarelativistic matter is negligible in comparison with 
the kinetic j2Hl and/or the potential [21] energy of the scalar field respectively. 

The inflationary solution can be obtained from the action if we impose restric- 
tions of the form 



which can be realized by an appropriate process of spontaneous symmetry breaking as 
described in the standard literature on this subject (see jSHj, jHOl IHl] and [H2])- The 
resulting model given below describes the inflationary expansion of an inhomogeneous 
universe in the limit g ^ oo (i.e. a oo), while the Kasner-like regime |Sj is obtained 
asymptotically approaching the singularity for — (i.e. a — oo). To derive such a 
solution we use the Einstein-Hamilton- Jacobi method since in the sense mentioned above 
this theory is the quasiclassical approximation to quantum gravity and also because it 
is computationally convenient. 

Let us consider the situation where in the potential (mass) term can be approx- 
imated as 



where A = A(x*) ~ const. These conditions impose peculiar restrictions on the degree of 
inhomogeneity of the scalar and gravitational fields and on the potential form of ^(0). 
Then the Einstein-Hamilton- Jacobi equations are 



If we separate the wave functional of the universe into its modulus and its phase, then the latter will 
satisfy the Hamilton- Jacobi equation containing an additional "quantum potential" . 




(14) 




(15) 




(16) 
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or explicitly 



E 



6a J 



+ n — 1) exp ina) A = 0. 



(17) 



The evolution of the residual variables remaining in the Kasner reference vectors 
1°" a are determined by the super-momentum constraint equations iJ^ = and can be 
expressed via the functions {Pr,P^). 

The solution of eq. (fT?)|l can be expressed in the form 



/(/3^a)= / S.KrP'^ + i-K +-ln 
Js [ nan 



where Ka{Kr,a) = i^/Xlr + ri{n — 1) A exp (na), K = ^/J2r are ar- 

bitrary "constant" functions of the spatial coordinates (i.e. independent of the time 
variable). Here S denotes the whole available spatial domain in which restrictions ()14|1 
are fulfilled. 

The signs ± before the square root correspond to the two possibilities for the variation 
of the local spatial volume in S (it depends on whether collapse or expansion of S is 
considered). Indeed, from p3p we find that the variation of a is determined by the 
Hamilton equation 

da ^ 2NP^ 

dt n {n — 1) exp {na/2) 

and N/ exp {na/2) > 0. In order to further simplify our analysis, we choose as a time 
variable the quantity a (i. e, -^a = 1), which implies the time gauge condition = 
n{n — 1) exp {na/2) / (— 2P„)) (since the lapse function should be positive by definition, 
we must have Pa < 0). 

Now according to the Hamilton-Jacobi method, we differentiate with respect to the 
quantities and then by putting the results equal to arbitrary "constant" functions, 
we find the solutions describing the trajectories of the system {-^ = (3q) 



(5'{a,x')= I3'^{x') + 



n\K\ 



In 



Ka + K 



(20) 



where (3^{x^) are arbitrary "constant" functions. In the asymptotic limit g ^ oo (i.e. 
a — s> oo =^ Ka — > c>o) the solution ()2()|1 transforms into the inflationary solution ob- 
tained in |27|, corresponding to the quasi-isotropization of the model since the parametric 
functions f3^ approach (exponentially) the "constant" values /?g(x*). In the opposite limit 
g ^ (i.e. a — oo =^ K) (j^UI) transforms into the generalized Kasner solution 

as it should, modified by the presence of the scalar field (see the |221) 



f3^{a,x') = (3l{x' 



— {a - ao) 



(21) 



where a^ renames the remaining constant terms. 
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We conclude our analysis by establishing the relation between our time variable a 
and the synchronous time T, which has a precise cosmological interpretation. These two 
variables are connected by the differential expression 

n (n — 1) exn (na/2) ^ , / x 

dt = Nda = ^ '—^ — —)da. (22) 

Now it is easy to see that from the Hamilton equation obtained by varying the action 
with respect to a we get, having fixed our time gauge (and remembering that H'^ = 0) 
the following asymptotic behaviors: 



a 



oo : Pq, oc —VXexp {na/2) , a — oo : Pa oc const < 0. (23) 



By substituting these relations into (f^ . we find the expected (familiar) asymptotic 
relations (which make evident the character of the two "opposite" regimes) 



a 



oo : a oc =^ ^ oc exp (^Ci\^t^ a — oo : a oc — lnC2t =^ y/g oc t 

(24) 

where Ci and C2 denote two constant values. 

The existence of this solution shows how the inflationary scenario can provide the 
necessary dynamical "bridge" between the fully anisotropic and the quasi-isotropic stage 
of the universe evolution. 



4 Misner-Chitre-like approach 

Though our solution is perfectly characterized by the above Misner-like variables, nev- 
ertheless to make precise the restrictions to be imposed on the initial conditions for the 
existence of such an interpolating regime, it is necessary to investigate a bit in detail the 
finite oscillating evolution to the singularity and therefore it is much more convenient 
to make use of the so-called Misner-Chitre-like variables. [51 In order to introduce these 
variables, the scale functions may instead be parameterized as follows (see Ref. [T!^ ) 

q'' = \nRl + MJng; ^M, = l, a = l,2,...,n (25) 

a 

where we distinguished a slowly varying function of time Ro, which characterizes the 
absolute value of amplitude of the metric functions [HI ESj and is specified by initial con- 
ditions (see below), from the anisotropy parameters Ma, which characterize the model's 
anisotropy; now the quantities Ing = — 2n In Rq and can be expressed in terms 

of the new set of Misner-Chitre-like variables r and {j = l,2...{n — 1)), as follows: 

\ng = -ne Ma[y^)=-[1 + — — ^ , 26 

I — y n \ 1 + J 

where A°'i is the matrix (fT^ . The Misner variables (3^ are related to ?/* by 

= -e--2yV (1 - y^) , a = \nRl - e"- (l + y^) / (l - y^) . 
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The parameterization ()2(j|l is defined within the domain — oo<r<oo,0<?/<l 



{y = y^jiy^y) and (with < (7 < 1) an appropriate choice of the function _Ro allows 
one to cover all of the classically allowed region of the configuration space using this 
parameterization. 

Within this choice of variables, the evolution of the scale functions is described by 
the action 

(27) 

where Pj/ and h denote respectively the conjugate momenta to y and r, while = 
^ (1 — 2/^)^ P| and the potential term 11 has the following structure 

a,b,c 

Here the coefficients Xabc (constructed by the spatial derivatives of the reference vectors 
£q) are slow functions of Ing, i.e. of the time variable, and characterize the initial 
intensity of the inhomogeneity field. When (7 -C 1, we can use the approximation of 
deep oscillations [H 121 u, in which the above potential is modeled by a set of potential 
walls, each of them having the form 

''^w ^ { • (29) 

As a result the whole potential becomes asymptotically (11 ^oo) independent of 
Kasner vectors, i.e. IIoo = J2^oo {<^a)- 

By solving the Hamiltonian constraint = in (jTTj) we define the ADM action 
pi] , reduced to the physical phase space, by a standard procedure which leads to the 
following (reduced) action 



'red 



/ (Py ■ % + Q^(l - HADM)d^xdT, (30) 



where 



Hadm = -h= + e-2-Q2 + n (31) 

is the ADM Hamiltonian and r now plays the role of the time variable (f = 1), once the 
gauge is fixed by Nadm = ^^H^^lf^ ■ For convenience we have also redefined the 
scalar field variables as follows q = (3"' and Q = Pn- 



^ It is a well-known result that, close enough to the singularity, the inhomogeneous Mixmaster evolution 
can be represented as a sequence of Kasner epochs (i.e. completely neglecting the potential term), 
while the transition from one epoch to the next one can be regarded as an instantaneous phenomenon 
(i.e. instant bounces against the potential walls). This well-established oscillating behavior is called 
the "deep oscillation approximation" and corresponds in our Hamiltonian picture to the replacement 
of the actual potential term with a set of infinite potential walls. 
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Thus, in the hmit r — oo, which corresponds to 11 ^ IIoo {y-') (the independence of 
the potential walls on the time variable r is the most relevant advantage in using Misner- 
Chitre-like variables), the system (jHUj) is nothing more than a point-like realization of 
a billiard on the (n — 1) —dimensional Lobachevsky space. By other words, in each 
point of the space, the system dynamics is isomorphic to the motion of a point particle 
on the (negative constant curvature) (n — l)-dimensional hypersurface formed by all 
the admissible values (the real gravitational degrees of freedom). On this domain the 
potential walls cut a region, which in dimensions n < 9 has a finite volum^ and therefore 
(in the absence of a scalar field) the resulting billiard exhibits strong mixing properties 
|3]. The role of scalar field (whose momentum, in this approximation, is simply an 
arbitrary "constant" function Q (x*)) in the evolution of metric functions consists of 
making geodesic lines on the billiard to be of a finite length and therefore of suppressing 
the chaotic regime (no other bounces against the potential walls take place). 

Thus, with this scheme in our hands, we are now able to easily specify the appro- 
priate inequalities characterizing different dynamical regimes. When approaching the 
singularity, the chaoticity can develop only in those regions of the universe where the 
energy of the scalar field is sufficiently small e~'^'^Q'^ ^ e^. We also observe that the 
condition for applying the approximation (j^^ corresponds to free motion in the allowed 
domain and therefore it can be written as follows 

> U. (32) 

From the condition that the approximation of deep oscillations (|^^ breaks at the 
moment (7 ~ 1 (r ~ 1), we find that the function Rq should be chosen as follows 
^2n-2 _ _^__g2T (^-^j^ere = E Aabd), so that in fact the inequality then just 
becomes g <^ 1 (r ^ 1). 

5 Dependence of the Universe Evolution on the Initial 
Conditions 

In this section we develop a synthesis of the interpolating regime constructed above, but 
in view of the different dynamical issues can take the universe in consequence of different 
initial conditions. All the considerations presented in the discussion below can be easily 
derived from the Hamiltonian function and equations associated with the reduced action 

As emphasized above in the presence of a scalar field, the final stage of the 
cosmological collapse ends with a monotonic Kasner-like behavior, i.e. the number of 
oscillations is always finite. From a phenomenological standpoint, if we consider an 
isotropic universe, this means that the effective "equation" of state for the gravitational 
waves (or the anisotropy parameters) is slightly softer then the one for the scalar field 

^ In the presence of dilaton-p-form fields billiards have finite volumes in all dimensions, e.g., see Refs. 
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(behaving like stiff matter, ctiaracterized hj e = p), i.e., e > p. To see ttiis explicitly, we 
may use some of the results derived above. 

Indeed near the singularity we have 11 — > IIoo [y) and the scalar field momentum does 
not depend explicitly on the time variable Q = const and therefore the same behavior 
characterizes the ADM "energy" of the anisotropy = e'^{y, Pg) = const (clearly from its 
expression, this quantity does not depend explicitly on r). However, the ADM energy of 
the scalar field depends on time = {QYe'^'^ and we get the relation tg/eh ~ e""^ ~ \ng. 
Thus in the limit r — oo {g 0) the scalar field dominates and the ADM Hamiltonian 
does not depend on the gravitational variables Hadm ^q, i-e. it turns out y = const, 
Pg = const (this corresponds to a stable Kasner regime approaching the singularity). 

On the other hand, during the expansion of the universe the role played by the 
scalar field through the evolution becomes smaller and smaller an, therefore the following 
various types of dynamical regimes can take place. 



5.1 The vacuum-type regime 

If the initial ADM energy of the scalar field is not very big, then the ADM energy of the 
anisotropy starts to dominate and this type of regime will be described by the oscillatory 
behavior with a frozen scalar field Q = const, q = const. The scalar field potential, which 
is in general always negligible near the singularity (i.e. at sufficiently high temperatures), 
in this case behaves like an effective cosmological constant V (q) = const. However if 
this constant remains small as compared to the ADM energy of the anisotropy during 
the whole period of applicability of the BKL picture, then this would result in the 
vacuum-type evolution. 

To define the moment r^n when the BKL approximation breaks down (i.e., g ~ 1), 
we consider the synchronous cosmological time, which is related to r by the equation 

n(n — 1) Rn\fq n , . 

dt = NADudr = -^—-L^^e-^^dT. (33) 

In the vacuum stage we have Hadm = = const, which gives us 

1 - \n{t/tin) 

where tin = "^^"^nel^" -Rq""^ — n(ri.-i)A2 rccall that during the expansion, Rq 

remains a constant in the oscillatory BKL regim, since as shown in Ref. [2111] A ~ Aoc"^. 

From the physical standpoint the moment Tin {din ~ 1) corresponds to the case where 
the characteristic scale of the inhomogeneity Lj„ matches the horizon size L/i, Lj ~ L/j, 
the so-called moment of origin of a stable background e.g., see |17j . 



5.2 The inflationary-type regime 

This kind of regime is realized when the scalar field potential grows until it is comparable 
with the ADM energy of anisotropy before the end of the BKL oscillatory regime. This 
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dynamical feature results in the so-called inflationary-like evolution, replacing the last 
Kasner epoch of the oscillatory regime (recall that if the scalar field potential remains 
small, then the evolution corresponds to the vacuum case). 

Under the assumption that the potential term gV [q) starts to dominate before the 
end of the Mixmaster r < ~ 1, then the ADM energy can be approximated by a 
function of the form 



I ADM 



^Jel + n{n-l) Rf'ge-'^^V (q) (35) 



(where we neglect the kinetic term of the scalar field and the spatial curvature term 
also). Thus the conditions for this regime to exit are expressed via the inequalities 

el « e2 ^ n (n - 1) Rl^g^e-'^^V (g) , (36) 

where Tc corresponds to the beginning of the inflationary evolution, i.e. the moment at 
which both terms at right hand side of (jHBj) have the same order, that is to say 

n{n-l)Rl-g,e-'^^V{q) 

and the inequality which expresses the applicability of the Mixmaster approximation 
reads gc -C gm ~ 1, or equivalently, in terms of the synchronous cosmological time, as 
follows from (jH^ 

tc < Un , (37) 

since tc is given by 

n — 1 , ^ 

^ V^- 

From a physical point of view the last inequality expresses the well known condition 
Li ^ Lc for the inflationary dynamics, where Lc is related to the field "Compton 
scale". We also recall that at later times, i.e., as t > tc, the inflationary evolution leads 
to the scalar curvature term in Hadm overtaking the kinetic energy of the anisotropy 
e\. However, we emphasize that both these terms will remain very small compared to 
the effective cosmological constant V{q). 

We note that the inflationary regime has a finite duration, due to a 'slow' evolution of 
the scalar field (the so-called slow-rolling phase), which is well described in the canonical 
literature on inflation [7j. The slow evolution phase of the scalar field can be found from 
the Hamiltonian equations 

TT^ = ^ ^ 0' (39) 

dr Hadm 

and 

d_ _ l n{n-l)ge-'^Rl-V{q)' 

Hence in the synchronous gauge (A^ = 1) dt/dr = "^""^-' -Rq y/g exp{—2T) / Hadm we get 

I, = '-9. (41) 
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= -R'^VaV (g)' , (42) 

where the initial conditions should be chosen so that 

Ql « ^y^' <n{n-l) Rl'^g^V (gj ■ (43) 

The rate of expansion is described by the equation dt/dr which can be rewritten as 
follows (we take \n.g = —ne~'^ {y = 0) and Hadm = \/n (n — 1) RQ"-ge~'^'^V (q)) 



d_ 
di 



V5-/^V^.V«^V5;exp(/*/^.). (44) 



5.3 The scalar-field-dominated-type regime 

Finally we consider another type of regime which takes place when the scalar field 
dominates during the entire Mixmaster approximation, i.e. ^ e^. In this case the 
solution has a monotonic behavior; wc can neglect the anisotropy functions from the 
very beginning and the evolution proceeds in the same way as in the case of isotropic 
models in the presence of a scalar field. 



5.4 Brief concluding remarks 

Wc conclude this section by emphasizing some relevant features of the interpolating 
solution we have obtained which help in getting physical insight into its cosmological 
setting. 

i) It is remarkable that the condition ensuring the inflationary scenario starts well 
inside the range of validity of the oscillatory regime is nothing more than the assumption 
at that moment Tc that the "Compton length" associated with the scalar field be much 
smaller than the inhomogeneous scale Lj„ of the universe; this condition is telling us that, 
as in the homogeneous case (when the condition for the existence of an inflationary stage 
reads in this same form), even in this general context the inflation can start only if the 
spatial gradients are sufficiently small. Thus this requirement for the validity of the BKL 
approach which is at the heart of our interpolating approach is just the inhomogeneous 
generalization of the well-know standard restriction for homogeneous inflation; this fact 
provides very strong physical support for the dynamical analysis developed here. 

However, we must mention that the condition so obtained is a point-like one and in 
principle, cannot take place in some spatial domain. It is worth noting that since wc 
require that inflation starts before the BKL regime ends, i.e. when the horizon size is 
still much less then the inhomogeneous scale {L^ -C L^), the dynamical regime derived 
above can (nevertheless) occur naturally in causal regions. 

ii) An important feature of the inflationary scenario consists of the frozen dynamics 
acquired by the anisotropy of the universe. In other words, during this stage of the 
evolution, the functions (3+ and (3- remain almost constant. As a consequence, the only 
evolving variable is now alpha, i.e. the volume of the universe; this fact makes it evident 
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that the dominant term during the inflation is Ae^" and any other term in the spatial 
curvature, although increasing like (at most) e^", becomes more and more negligible. 

iii) Finally we emphasize that the results obtained above do not depend on the 
particular form taken by the scalar field potential, as long as it realizes a spontaneous 
symmetry breaking process and a sufficiently long phase of slow- rolling 

It is also worth noting that the compatibility of the inflationary scenario with our 
actual cosmic phenomenology seems to be confirmed by recent observations on the mi- 
crowaves background radiation thus allowing, on the basis of the analysis developed 
here, very general dynamical behavior for the primordial phases of the universe evolution. 

The referee is thanked for valuable contributions in improving the form of this paper. 
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